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In this paper the dynamic characteristics of spiral-grooved seals are 
theoretically obtained by using the Navler-Stokes equation. First, with the 
Inertia term of the fluid considered, the flow and pressure In the steady state 
are obtained for the directions parallel to and perpendicular to the groove. 
Next, the dynamic character Is obtained by analyzing the steady state by ana- 
lyzing the labyrinth seal. 

As a result, the following conclusions were drawn: 

(1) As the land width becomes shorter or the helix angle decreases, the 
cross-coupling stiffness, direct and cross-coupling damping, and add mass 
coefficients decrease. 

(2) As the axial Reynolds number Increases, the stiffness and damping 
coefficients Increase. But the add mass coefficient Is not Influenced by the 
axial Reynolds number. 

(3) The rotational Reynolds number Influences greatly the direct and 
cross-coupling stiffness and direct damping coefficients. 

(4) As the journal rotating frequency Increases, the leakage flow 
decreases. Therefore zero net leakage flow Is possible at a particular 
rotating frequency. 


INTRODUCTION 

High-performance pumps, ( 1 . e . , those operating at high rotating speed and 
high pressure) are used In chemical plants, rocket engines, etc.. These pumps 
sometimes yield nonsynchronous vibration that Is Induced by the journal bear- 
ings or noncontacting seals. The Instability of rotors supported by a journal 
bearing has been studied very well. However the noncontacting seal, which 
sometimes Induces a nonsynchronous vibration or Instability to the pumps, has 
not been Investigated. 
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The dynamic characteristics of some noncontacting pump seals (the annular 
plain seal, the tapered seal, and the stepped seal) have been theoretically 
investigated by Black (refs. 1 and 2), Childs (refs. 3 and 4), and Yang et al . 
(ref. 5), but there have been no investigations using the theory of fluid 
dynamics on the parallel-grooved and spiral-grooved seals. Childs studied the 
parallel-grooved seal by using only a rough approximation. On the other hand, 
the parallel-grooved labyrinth seal has been investigated by Kostyuk (ref. 6) 
and Iwatsubo, et al. (ref. 7). For the investigation of the spiral-grooved 
seal, we studied steady-state characteristics from the viewpoint of leakage. 
This study was based on investigations of the spiral journal bearing; that is, 
Whipple (ref. 8) presents a basic idea to analyze the characteristics of the 
thrust spiral-grooved bearing, Vohr (refs. 9 and 10) and Passera (ref. 11) 
present an approximate method using creeping flow analysis (which assumes that 
the groove number is infinitely large), and Zuk (ref. 12) analyzes the static 
characteristics of the spiral-grooved seal by solving the Navier-Stokes equa- 
tion with a finite difference method. But all these analyses are for the 
static characteristics and there are no investigations of the dynamic 
characteri sties . 

This paper presents an analytical method to obtain the dynamic character- 
istics by solving a Navier-Stokes equation with the perturbation method. First 
the steady-state flow and pressure distribution in the axial direction are 
obtained by considering the pumping effect due to the spiral groove. Then the 
dynamic characteristics are obtained at the steady-state condition, and the 
dynamic force is represented by the matrix form. 


SYMBOLS 

C iq mean clearance 

D seal diameter 

F flow Induced force 

H thickness of fluid film 

I thread number 

s 

L seal length 

L groove width 

9 

land width 
L L 7 /L 

zg z g 

L $ number of lands for one thread 

NN number of land and groove sections 

P pressure 

AP pressure difference 

Q leakage flow rate 
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R journal radius 

R axial Reynolds number 

d 

R g Reynolds number 

R r circumferential Reynolds number 

T depth of ditch 

t time 

u circumferential fluid velocity 

w axial fluid velocity 

x,y,z coordinates (see fig. 2) 
a spiral angle 

c perturbation coefficient 

e divergent flow angle 

\ friction coefficient 

v viscous coefficient 

}■ loss factor 

£,n,C coordinates (see fig. 2) 
p density 

t shear stress 

Subscripts : 

a axial direction 

c radial direction 

d ditch (mainly used for vortex in ditch) 

ex exit 

f between clearance flow and cavity flow 

g groove 

in inlet 

j journal 

l land 
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m mean velocity 

r circumferential direction 

s casing 

(”) time mean 


GOVERNING EQUATION AND MODELING OF SPIRAL-GROOVED SEAL 
Governing Equation 

Figure 1 Illustrates the geometry of the spiral-grooved seal. Under the 
usual assumptions for problems of through flow across annuli with fine clear- 
ances, the continuity and momentum equations are represented as follows 
(ref. 13): 

Continuity equation, 

_3H + 3 (£V + 3 (HWjjj) _q ( . 

3 t 3a: 3 s 


where u^ and w m are mean fluid velocity components In the tangential 
and axial directions. 

Momentum equation In the x-dlrectlon, 


p H{ lUE + u m ^HL + }= -H— + T* r 


3 1 


3 w 


3s 


3a: 


S 

0 


( 2 ) 


Momentum equation In the z-dlrectlon, 


pff{ l + }= -H— + 


3 1 


dx 


3s 


3 z 


T ta 


E 

0 


(3) 


Modeling 

There are three kinds of spiral-grooved seals, (1) those with the groove 
on the journal (fig. 1), (2) those with the groove on the casing, and (3) those 
with grooves on both journal and casing. In this paper the seal with the 
groove on the journal Is analyzed. Figure 2 Is an expanded figure of the 
spiral-grooved seal. Configuration parameters are Indicated on this figure; 
these are spiral angle a , land width L ^, groove width Lg, seal diameter 
D, and thread number Is- These parameters are related as follows: 
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( 4 ) 


where L s Is number of lands for one thread In the axial direction and 
L^g = L^/Lg . For the analysis two coordinates are used: one Is the 

n-y-f coordinate, which Is used for the static analysis and the other Is the 
x-y-z coordinate, which Is used for the analysis of dynamic characteristics. 

For the analysis the following are assumed: 

(1) The fluid Is liquid and noncompresslble. 

(2) Flow In the land In the n-dlrectlon Is assumed to be a flow between 
two parallel plates, and flow In the groove Is assumed to be a flow In a rec- 
tangular cross section and Is approximated to the flow In a circular tube. 

(3) Flow In the C-dlrectlon diverges with the angle e and goes to 
the next land. 

(4) The vortex Is formed In the ^-direction of the groove. But the 
heat energy is negligibly small. 

(5) When the journal deviates a little from the center, the streamline In 
the groove deviates In the same manner. 


DERIVATION OF SHEAR STRESS IN MOMENTUM EQUATION 

The shear stress term of equation (3) Is derived for the spiral-grooved 
seal. The friction coefficient of the annular seal Is represented by a 
Reynolds number (ref. 14). But It Is very difficult to represent the friction 
coefficient of the spiral-grooved seal In the same way because the groove Is 
spiral. So the friction coefficient equation derived by Hlrs (ref. 15) Is used 
in this analysis. The equation of the friction coefficient between two plates 
Is represented as 


- 0.2 5 

X = 0 . 066 Re 


(5) 


where the velocity used In the Reynolds number Is the equivalent mean velocity 
that Includes the pumping action of the spiral-grooved seal. This equivalent 
mean velocity Is used to obtain the shear stress. 


Land 

The shear stresses of the casing and journal parts In the axial direction 
( z-di rectlon) T £sa* T £ja are represented by the formula for the flow 
between two parallel plates, 
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T lsa X lja 


— pX wl 
2 la Im 
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Therefore the shear stress term of equation (3) becomes as 

I H 

= -2 


' ta 


X lja X lsa 


2 P ^la W lm 


(7) 


Groove 


In the groove, flow Is divided Into two regions: the jet flow region and 

the vortex region. The shear stress of the casing part of the jet region In 
the axial direction Tg Sa is represented by the formula of the flow 
between two parallel plates, 


T 

gsa 



ga gm 


( w h r°- 2s 

X = 0.066 -i £ 

ga \ v J 


( 8 ) 


It Is assumed that the energy loss In the groove vortex region Is represented 
by the friction loss between the clearance flow and vortex flow. Assuming that 
the friction between the vortex and the jet flow \f Is 0.1 and that the 
axial velocity component of the vortex w,j m Is set to one-half of the jet 
flow w gm (ref. 16), the shear stress of the flow Tgf a becomes 

Va - - t P V V - "a, ! ‘ ■ < 9 > 


Therefore the shear stress term In the groove Is obtained from equations (8) 
and (9) as 


T 




T 

gsa 


- -j p 0.25X 


f gm 


— — pX w 
2 ga 


2 

gm 



+0.25X 



( 10 ) 


DERIVATION OF STEADY FLOW AND STATIC CHARACTERISTICS 
Axial and Circumferential Velocity of Steady Flow 

As described before, the n-y-C coordinate system Is used for the ana- 
lysis of static characteristics. References 10 and 11 describe investigations 
of the static characteristics of the spiral-grooved seal. In these Investiga- 
tions, the residual seal pressure Is obtained when seal leakage becomes zero 
because of the pressure Induced by spiral pumping action. This residual pres- 
sure Is a function of the seal configuration, and the seal coefficient can be 
represented by the nondlmenslonal form. Boon et al . (ref. 17) obtained the 
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seal coefficient of laminar flow, and Mori, et al. obtained the same by exper- 
iment. The seal coefficient S-C-iamlnar for the laminar flow Is 


S.C. 


laminar 




k 3 ( 1 + tan 2 a)+L . ( 1-L , ) (k 3 -1) 2 tan 2 a. 
l <3 ig 

L - (1-L. ) (k 3 -l) ( k- 1) tana 
ig ig 


(ii) 


where k = Hg/H e . THe pressure Induced by the spiral pumping action 
APlamlnar for the laminar flow Is 


AP 


laminar 


6\iHuL L . (1-L. ) (< 3 -l) (k- 1) tana 

l£ Ig 

Cl k 3 ( 1 + tan 2 a.) +L . (1-L. ) (< 3 -1 ) 2 tan 2 a 
lo Ig Ig 


( 12 ) 


Vohr (ref. 10) compared the seal coefficient for turbulent flow with that for 
laminar flow by experiment. He obtained the pressure ratio of the seal 
coefficient of the turbulent and laminar flow as 


C - — turbulent _ 0m Qj 5g % 0.7 7 8 

* AP r 

laminar 


Then the pressure AP^ ur b U i e nt Induced by the spiral pumping action In tur- 
bulent flow Is 


0. 09Sl)iRuLR 


1.7 7 8 


AP 


turbulent 


C lo 

L. (1-L. ) (< 3 -l) (k-1) tana. 

lg_ Ig 

ic 3 (1 + tan 2 a)+ L. (1-L. ) (k 3 -1) 2 tan.a 
tg Ig 


(14) 


Since the pressure difference equation (14) acts to resist the pressure differ- 
ence In the seal AP, the pressure difference In the seal becomes smaller than 
the pressure difference between both seal ends. This pressure difference Is 
called the apparent pressure difference AP'; that Is, 


AP" - AP - AP 


turbulent 


(15) 
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This apparent pressure difference acts on the seal and on the fluid flow In the 
seal. This flow Is separated Into n- and c-dlrectlons and the fluid veloc- 
ities in each direction and In each part are derived In order to calculate the 
dynamlc-f low-induced force In each stage. Then total force Is obtained by sum- 
ming up the lateral force of each stage. 

Fluid velocity parallel to groove direction for land . - The relation 
between the pressure and velocity within the two parallel plates Is also con- 
sidered In the land region; 


A P' 


1_ 

2 


P ( l + £> 


lr\in 


)w 


2 

lr[o 




2 

It)o 


2L 


C, sina 
lo 


+ 2 (16) 

The first term on the right-hand side of equation (16) represents the Inlet 
pressure loss. The Inlet pressure loss of the land Is larger than that of the 
groove because the Inlet clearance of the land Is smaller than that of the 
groove and because the fluid flows Into the groove. The second term on the 
right-hand side of equation (16) represents the pressure loss due to wall 
friction. The third term represents the outlet pressure loss. The value X^ n 
Is obtained using equation (5), and 5j n 1n Is the Inlet pressure loss 
coefficient, which Is 0.5 for the first stage and Is represented as follows 
after the second stage: 


£i n in = 1 +0. 8246 2 - ( 1+0.8246 i )(H 2 /H 1 

61 = 1 -95 ( w 10 H-| / v )-°- 43 

6 2 = 1 -95 ( w 10 H 2 / v r 0,43 

Therefore If equation (14) Is Iteratively calculated so that w^ n o 
function X coincides with w^ n o In equation (16), the flow velocity of the 
land w £ n Is obtained. 

Fluid velocity parallel to the groove direction for groove . - Since the 
groove Is deep enough, the groove Is considered to be a rectangular pipe and 
can be approximated as a circular pipe. The equivalent radius of the rectan- 
gular pipe MG (ref. 18) Is 



MG 


L ( C 7 +T) 
g__Lo 

2 (L +C, +T) 
g lo 


( 18 ) 


The friction coefficient of the turbulent circular pipe Xj (ref. 18) Is 
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Therefore the apparent pressure difference AP' 


is written as follows: 


A P' = 


2 g n in gno 


+ — A-.— 

2 d gno sma MG 



1- 


C , +T 
lo 

C, + 

lo 


w‘ 


gno 


( 20 ) 


In this equation the first term on the right-hand side represents the Inlet 
pressure loss of the groove. This Inlet pressure loss Is smaller than that of 
the land, because the groove has a larger cross section than the land and 
because the liquid flows mainly Into the groove. The second term represents 
the pressure loss due to wall friction, and the third term represents the out- 
let pressure loss of the seal. After Sg n i n Is obtained by using equa- 
tion (17), the velocity In the groove In the n-dlrectlon wgr,o Is similarly 
obtained by Iteratively calculating equation (20). 

Fluid velocity perpendicular to groove direction . - The flow In the 
C-dlrectlon Is derived by a method similar to that used to obtain the steady 
flow of the parallel-grooved seal. Conventional evaluation of the loss of the 
groove region was not clear. Yamada (ref. 19) obtained the friction coeffi- 
cient of both the land and the groove by experimental methods. According to 
his results, as the groove region Increases In relation to the land region, the 
friction coefficient Increases as shown In figure 3. This friction coefficient 
Is represented by using the results of reference 19 

-0.24 7 -L , 

X = 0.26'Re • 3. 31 9 (21) 

From the pressure drop relation the following Is obtained: 


\ ( [/ i+L } — X-L- X L 

i g ii g g 


( 22 ) 


where \i Is a friction 
the friction coefficient 
and (22), 


coefficient of 
of the groove 


flow between two parallel plates. Then 
\ a Is obtained from equations (21) 


X = 0.26 Re 

9 


-0.24 


l-l 7 L 7 
3. 31 19 +~r^( 

9 


1-L - , 

3.3i 19 -1 n 
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Therefore the apparent pressure difference In the f-dlrectlon is obtained 
by summing each stage of the land and groove, respectively, 


LP' 


1 % 2 7 ou 
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9 
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In this equation the first term on the right-hand side represents the Inlet 
pressure loss of the seal, the second and the third terms represent the pres- 
sure loss due to wall friction, and the fourth term represents the outlet 
pressure loss of the seal. The term In this equation Is defined In 

equation (5), 5£ln Is defined In equation (17) and NN Is the number 
of the land and groove sections In the ^-direction. Then the velocity In 
the c-dlrectlon w ^g i S obtained by a similar Iterative calculation. 


Fluid velocity In circumferential and axial directions . - The steady-flow 
velocity In the x- and z-dlrectlons Is obtained from the previously calcu- 
lated steady-flow velocities w ^ n0 , w gil0 , and w^^ 0 by translating the coordi- 
nate system as 



where u '^ 0 and u'g 0 are the velocities relative to the journal. There- 
fore the absolute velocities In the circumferential direction become 


u 


lo 





(26) 


Leakage Flow Rate 

Because the leakage flow rates of the land and groove are different, they 
are considered separately. The leakage flow rate for the land Is represented 
by the relation of the cross-sectional area and velocity, 


w 7 2-nrdr - ^C. I 2 (R+T)+C. w . (27) 

lo lo\. to) Ig lo ' ' 


The leakage flow rate for the groove Is derived by separating the clear- 
ance flow from the cavity flow. The leakage flow rate of the clearance flow 
In the axial direction Is represented as 


«l ' r s * I 


L _y\ K+M lo 


s ,/ R+T 
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The leakage flow rate in the cavity Is represented by the axial component of 
cavity flow as 


1-L R+T 

Q . = I x — - w _ _2nrdr 


' 9 d 


s J R 


go 


« Tr (T+C- )(2RtT+C 1 )(1-L- )w 

lo to Ig go 


(29) 


Therefore the total leakage flow rate for the groove Is 


Q 


g 




(30) 


and the total leakage flow rate Q becomes 


Q = Q, + Q 
4 g 


lo 

■hi 




- 1,C ' j2(!1+I)+C lo} L l g u lo 


C* L 


7 i 

~L ( 2C , +L tanB)w sina-h -^2. w 
2 g lo g gqo tan 8 lZo 
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lo lo Ig go 


(31) 


Figure 4 shows the leakage flow rate for three types of seal and for 
L/D =1.0, a rotating speed of 4000 rpm, and a radial clearance of 0.5 mm. The 
leakage flow rate for each seal type Increased as the pressure difference 
Increased. The leakage flow rate for the parallel-grooved seal was less than 
that of the land seal. For the spiral-grooved seal. If the spiral angle was 
small, the leakage flow rate was less than that of the parallel-grooved seal, 
but If the spiral angle was large, the leakage flow rate was greater than It 
was for the parallel-grooved seal. For the spiral-grooved seal, If the spiral 
angle was small and the land width was large, or if the spiral angle was large 
and the land width was large, the leakage flow rate decreased. For this reason 
the groove is Important for the screw pumping action. 


165 


Figure 5 compares the calculated leakage flow rate with experimental 
results (ref. 20). The groove was on both the journal and casing In this 
experiment. In the low-speed range the calculated and experimental values were 
very close, and In the high-speed range the leakage flow rate tended to zero 
because of gas Ingestion phenomenon. 


DERIVATION OF PRESSURE DISTRIBUTION 


To linearize the equations, a first-order approximation Is performed. It 
Is assumed In the linearization that the center of the journal coincides with 
the center of the casing and that the journal perturbs close to the center. 
Then the values of fluid film thickness H, pressure P, mean velocity In the 
z-dlrectlon w m , and mean velocity In the x-dlrectlon u m are written as 
follows: 


H = C + c\ b 
o 


w = w + zw , 
mo 1 


P = P + 

o 


u 


m 


u 

o 



(32) 


where C 0 , P 0 , w 0 , and u 0 are the steady-state values and y, P-j , W] , and 
u-j are the perturbations. The perturbation term In the circumferential 
direction Is neglected because the spiral angle a Is usually small and u-| 
Is sufficiently smaller than W] . 


Pressure Distribution In Steady-State 

The steady-state pressure distribution Is obtained by substituting equa- 
tions (7), (10), and (32) Into equation (3) and taking the zeroth-order approx- 
imation. The pressure gradient of the land In the axial direction Is given by 


3P 


lo 


3 3 


0.132 


lo 


' 2 QW lo 


M 2 ] 


0.2 S 


(33) 


The pressure gradient of the groove In the axial direction Is given by 


3P 
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where Cg 0 and w qo are functions of z and are represented as 


y go 


C = C , + ztanBcosa 

go lo 
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C lo W lo 


g° c lo +*tanBcos* 


(35) 


Therefore equation (31) Is written as 
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( 36 ) 


3P 



1 ^ Iq W Iq) Xsg 

2(C~ +ztanQcosa)3 

to 


pC - u V- tanBcosa 
Lo g_o_ l_o 

( C i o +ztanQcosct) 2 


where \ sa Is the apparent friction coefficient because Its dimension Is the 
same as friction; that Is, 


X = X + 0.25X-, ~2tanBoos a 
sa ga f 


(37) 


Dynamic Pressure Distribution 

The dynamic pressure distribution Is obtained by substituting equa- 
tions (7), (10), and (32) Into the equation of momentum (3) and taking the 
first-order approximation. For the continuity equation (1) the first-order 
approximation Is also obtained In the same manner. The momentum equation In 
the axial direction for the land Is 


3P 
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i|> 3P , /C, w 7 \ w -‘* 
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(38) 


The continuity equation Is 


3u) 


U = , M 

3 z C^ o \%t to 3 xj 


(39) 


Pressure distribution Is obtained by solving equations (38) and (39) simulta- 
neously. Equation (39) Is the first-order ordinary equation In w^i , and It 
Is solved by setting the boundary condition w^i = w^i(0) at z = 0. 
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11 


= w 


11 ( 0 ) 


to 




(40) 


where 




11 J1 

dt to 3a 
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Substituting equation (40) into equation (38) and setting the boundary condi- 
tion P £•) = P^i ( 0) at z = 0 yield a solution of equation (38). In this 
boundary condition P£i(0) is the perturbed pressure at the land inlet which 
is obtained by a first-order approximation of the equation 

AP = - (1/2)pC-|W m 2 (41) 


that is. 


P ll(0) 


P ^l U lo W ll(0) 


(42) 


The term w^-|(0) is obtained by neglecting the perturbed pressure at the 
outlet of the land, (i.e., z = L^/cos a, P^i (0) = 0). Then the dynamic 
pressure of the land P-p is obtained as 



2C * 
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where SPL1 and SPL2 are as written in appendix A. 

The momentum equation in the axial direction for the groove is 
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The continuity equation is 
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Substituting equation (36) Into equation (45) and setting the boundary condi- 
tion Wg-j = Wg i ( 0 ) at z * 0 yield the solution of equation (45): 


w 
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to 
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w i o 
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where 


f (Xjt) 
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Substituting equations (36) and (46) Into equation (44) and setting the bound- 
ary condition Pg-j = Pg-| ( 0) at z = 0 yield the solution of equation (44), 
where Pgi(O) Is the pressure loss at the Inlet of the groove: that Is, the 
outlet loss of the land. The outlet loss of the land Is represented by 


LF lex 'i P ^ U lo 

where 1s ‘ t ^ e outlet loss coefficient 

S 2 - ('•Cjrtjfne) 

The first approximation of equation (47) Is represented as 

p g l( Zn ) = -PC 2 w 10 w gl( z n> 


(47) 


(48) 


(49) 


The term Wg-|(0) Is obtained by neglecting the perturbed pressure at the 
groove outlet; that Is, z = Lg/cos a, Pgi = 0. Then the dynamic pressure 
of the groove Is obtained as 
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where SPG1 to SPG9 are represented In appendix B. 


DYNAMIC FORCE 

The small displacements In the X- and Y-dlrectlons at the center, AX 
and AY respectively, are as follows: 

eip = - AX cos <j> - AY szn<t> (51) 

The flow-induced forces In the X- and Y-dlrectlons at the beginning of a land 
4> n , at the beginning of a groove 4> n +l/2* an ^ at *be beginning of the 
next land section <j> n +l are represented by equations (52) and (53). 
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This force may be represented In the following matrix form: 


K 

K 


x 


c 
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X 


M „ 
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XX 

XY 



+ 

XX 

XY 



•f 

XX 

XY 



k yx 

k yy 

4 


Y 


C YX 

C Y Y 


Y 


m yx 

w yy 


Y 


(54) 


The following coefficients are obtained from equations (52) to (54). The 
dynamic coefficients are different for the land and the groove. 
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where ISPL2 to ISPG9 are represented In appendix C. 

The coefficient matrix represents the summation of the coefficients of the 
land and the groove. 


NUMERICAL EXAMPLE 

Data for the seal model are shown In table I. Figure 6 shows the relation 
between the axial Reynolds number (pressure difference between the Inlet and 
outlet of the seal) and the spiral angle a, where L/D = 1.0, rotating speed 
Is 3000 rpm, and pressure difference between Inlet and outlet Is 0.49 to 
4.9 MPa. As Reynolds number In the axial direction Increases, K xx , K yx , C xx , 
and C yx Increase, but M xx Is almost constant. For the 3-thread 2.6° 
spiral-angle spiral-grooved seal and the 20-thread 17.66° spiral-angle spiral- 
grooved seal, the dynamic coefficient decreases as the spiral angle decreases. 

Figure 7 shows the effects of the Reynolds number In the axial direction 
and the ratio of the land width to the groove width. It Is known that as the 
land width Increases, the coefficients K xx , K yx , C xx , C yx , and M xx Increase 
and that. If the land Is narrow, the spring coefficient Kxx becomes 
negative. 

Figure 8 shows the effect of the circumferential Reynolds number for two 
kinds of spiral-grooved and parallel-grooved seals, where L/D = 1.0, the pres- 
sure difference between the Inlet and outlet Is 0.49 to 4.9 MPa, and rotating 
speed Is 2000 to 8000 rpm. Coefficient K xx of the spiral-grooved seal shows 

a negative value for a low circumferential Reynolds number, but that of the 
parallel-grooved seal does not show a negative value until a high circumferen- 
tial Reynolds number Is obtained. The reason Is that the pressure difference 
between the Inlet and outlet of the spiral groove seal affects the circumferen- 
tial velocity; that Is, the pressure difference and the high rotating speed 
Induce the pumping action, and the apparent pressure difference becomes small. 
The characteristics of the coefficients K xy and C xx can be Illustrated In 
the same way. The values of the coefficients K yx , C xx . C yx , and M xx for the 
spiral-grooved seal are larger than those for the parallel-grooved seal. 


CONCLUSIONS 

From this study of the dynamic characteristics of noncontacting spiral- 
grooved seals, the following conclusions are drawn: 

1. Except for K xx , coefficients become small as the spiral angle and 
the land width decrease. 

2. As the axial Reynolds number Increases, coefficients K xx , K yx , C xx , 
and Cy X become large but M xx remains constant. 

3. As the circumferential Reynolds number Increases, K xx decreases, K yx 
and C xx Increase once and then decrease, C yx Increases, and M X x Is 
almost constant. 
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4. The coefficient K xx becomes negative when the circumferential 
Reynolds number Is large and the axial Reynolds number Is small. 

5. Leakage flow decreases as the rotating speed becomes large. If the 
spiral angle Is small, leakage flow decreases as the groove width becomes 
large, and If the spiral angle Is large, the leakage flow decreases as the 
groove width becomes small. 

6. Comparing the coefficients of the spiral-grooved seal with those of the 
parallel-grooved seal, yields the following conclusions: 

a. K xx of the spiral-grooved seal Is smaller than that of the 
parallel-grooved seal. 

b. K yx and C xx of the spiral-grooved seal are smaller (larger) 

than they are for the parallel-grooved seal for low (high) circumferential 
Reynolds numbers, respectively. 

c. C yx and M xx of the spiral-grooved seal are larger than they 
are for the parallel-grooved seal. 

d. Leakage flow of the spiral-grooved seal Is larger (smaller) than 
that of the parallel-grooved seal for low (high) rotating speeds, respectively. 


APPENDIX A 


«i 




2 C 


lo 


* W lo° 





cos a 


‘’ll " 0U la K l 


‘>11 ‘ IcJ^lc 0 - 01 ^— ) 


SPL1 = z 

JLln - n z) 

h 

‘ q^ q H q l2 J 

cos a 

SPL2 = 2 2 

-T^u - q l 2-*} 

( h 1 
\cosa / 


174 


APPENDIX B 
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TABLE I. - DATA FOR NUMERICAL CALCULATION 


Working fluid water 

Fluid temperature, K 293.15 

Density, p, Kg/m 3 9.982xl0 2 

Viscosity, y, mPa s 1.009 

Kinematic viscosity, v, m2/s 1.006x10-6 

Journal radius, R, mm 100.0 

Seal radial clearance, C^ 0 , mm 0.4 

Seal length, L, mm 200.0 

Groove depth, T, ran 1.18 

Helix angle, a, deg 1.30 to 17.66 

Journal rotating frequency, N, rpm 2000 to 8000 

Pressure difference, AP, MPa 0.49 to 4.9 


p in Y////////////A r, 



W/////////ZA 

Figure 1. - Spl ral -grooved seal. 
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Figure 2. - Illustration of spiral-grooved seal. 



Figure 3. - Friction coefficients. 




Figure 4. - Leakage flow versus pressure difference 



ROTATIONAL FREQUENCY N (rpm) 


Figure 5. - Leakage flow versus rotating speed. 
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(c) Damping coefficient, C xx . 
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(d) Damping coefficient, C . 

yx 

Figure 6. - Continued. 
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Figure 7. - Concluded. 
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STIFFNESS COEFFICIENT -Kyx (N/m) STIFFNESS COEFFICIENT Kx 



10 

ROTATIONAL REYNOLDS NUNBER Rr 

(a) Stiffness coefficient, K . 

' xx 
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ROTATIONAL REYNOLDS NUNBER Rr 

(b) Stiffness coefficient, K . 

Figure 8. - Effects of circumferential Reynolds number 
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ROTRTIOhfiL REYNOLDS NU1BER Rr 

(d) Damping coefficient, C . 

yx 

Figure 8. - Continued. 
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